module my_array
   use kind_parameter
   implicit none
   private
   public:: f

contains

   subroutine f(r)
      real(dp), intent(out) :: r(:)
      integer :: n, i
      n = size(r)
      do i = 1, n
         r(i) = 1.0_dp / i**2
      end do
   end subroutine f

end module my_array
